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We study the fidelity susceptibility in an antiferromagnetic spin-1 XXZ chain numerically. By 
using the density-matrix renormalization group method, the effects of the alternating single-site 
anisotropy D on fidelity susceptibility are investigated. Its relation with the quantum phase tran¬ 
sition is analyzed. It is found that the quantum phase transition from the Haldane spin liquid to 
periodic Neel spin solid can be well characterized by the fidelity. Finite size scaling of fidelity sus¬ 
ceptibility shows a power-law divergence at criticality, which indicates the quantum phase transition 
is of second order. The results are confirmed by the second derivative of the ground-state energy. 

We also study the relationship between the entanglement entropy, the Schmidt gap and quantum 
phase transitions. Conclusions drawn from these quantum information observables agree well with 
each other. 

PACS numbers: 03.67.-a,05.30.Jp 


I. INTRODUCTION 

In the last two decades, quantum phase transitions 
(QPTs) in quantum spin chains have attracted con¬ 
siderable interest both in experimental and theoretical 
research[T]. Among them, the spin S' = 1 antiferromag¬ 
netic Heisenberg chain has been extensively studied by 
many authors[^-[l], whose ground state is termed the Hal¬ 
dane phase. Such phase has a peculiar nonlocal string 
order in which the spins with |Sz = ±1) are arranged an- 
tiferromagnetically if the sites with jSz =0) are skipped. 
As we know, the QPT has been traditionally described 
based on the behavior of expectation values of local op¬ 
erators and two-point correlators in accordance to the 
standard Ginzburg-Landau theory. Unlike the case of 
spontaneous symmetry breaking, the absence of the long- 
range order in Haldane phase is accompanied by a hidden 
Z 2 X Z 2 symmetry breaking. Furthermore, the Haldane 
state is gapped between a spin-singlet ground state and 
a spin-triplet excited state, indicating spin-1 Heisenberg 
model is sharply different from spin-1/2 counterpart. 

Experimentally, a few quasi-one-dimensional Hal¬ 
dane chain compounds have been investigated, such as 
Y2BaNi05[S],CsNiCl3[BHH]) Dy2BaNi05[5], Nd2BaNi05 
[laiii] were considered to realize 5" = 1 Haldane systems 
with a magnetic gap in the excitation spectra. However, 
an ideal one-dimensional (ID) spin-1 system is rare in 
real materials, usually followed by the interchain inter¬ 
actions and magnetic anisotropy, which may partially or 
completely suppress the excitation gap and thus lead to 
an observation of long-range order in a quantum disor¬ 
dered magnet. The strength of single-ion anisotropy D 
was retrieved from inelastic neutron-scattering on quasi- 
one-dimensional spin-1 chain compound ANi2V208 (A 
=Pb and Sr) [T^] , electron spin resonance (ESR) study on 


PbNi->VQOg|13j. multifrequency ESR transmission spec¬ 
troscopy on single crystals SrNi2V208[Ill[Ii]. With large 
single-site anisotropy (D>0), the Haldane ground state 
changes to the large-D state without explicit order, where 
all spins are confined to the configurations |S'z = 0). 
On the other hand, the reverse single-site anisotropy D 
{D<0) also drives the Haldane state into the Neel state. 
Note that in recent years oxide heterostructure appears 
to be an appealing discovery platform for emergent elec¬ 
tronic states[TH]. More explicitly, by suitable mechan¬ 
ical, electrical, or optical boundary conditions at oxide 
interfaces, intersite couplings and novel magnetic states 
can be manipulated. For instance, a possible scenario 
for realization of ID zigzag chain in layered structures 
of transitional metal oxide was proposed |17j. The re¬ 
duced symmetry at the interface strongly modifies the 
magnetocrystalline anisotropy at the interface. As a con¬ 
sequence, amplitudes and even the signs of the effective 
single-ion anisotropy can fluctuate. Hence the large in¬ 
terest arises when considering the anisotropy effects, that 
is, the exchange anisotropy and the easy-plane staggered 
single-ion anisotropy coexist in spin-1 chains. 

Recently, various exogenous approaches inherited from 
quantum information to understanding quantum many- 
body systems exploit the curvature measures of the 
ground states. Much effort focused on the study of QPTs 
in spin chains in terms of entanglement entropy (EE) [18j . 
Another concept referred to simply as fidelity suscepti¬ 
bility (ES), which measures the changing rate between 
two states, will diverge at the critical points[T9]. The 
ground-state EE and the FS were deemed plausible to 
qualify QPTs in strong correlated systems [201150] . It is 
due to QPTs are intuitively accompanied by an abrupt 
change in the structure of the ground-state wave func¬ 
tion. This primary observation motivates researchers to 
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use the EE and the FS to predict QPTs. The scaling rela¬ 
tion of the FS and the EE proposed recently is verified for 
the spin-1 XXZ spin chain with a single-site anisotropy 
term|31j. Through a proper finite-size scaling analysis, 
the results from both the FS and the EE, agree with the 
findings in the previous results (32] • The effect of spatial 
modulation of site-dependent anisotropy Di in the S = 1 
Heisenberg chain was studied using perturbation theory 
and exact diagonalization on small-size system [HH]. Thus, 
in order to figure out the QPT in ID spin-1 Heisenberg 
chains with exchange anisotropy and alternating single¬ 
site anisotropy, it is better to consider these two quantum 
information observables again. 

In the present paper, we calculate the ground-state FS, 
the EE and the Schmidt gap in ID spin-1 XXZ chains 
with alternating single-site anisotropy, and make them 
ideal tools for searching the phase transition points. The 
paper is organized as following. In Sec an introduc¬ 
tion is provided. The Hamiltonian is shown in Sec [H] 
The measurements and details of method to obtain the 
ground-state properties are introduced in Sec. |HI| In 
Sec |IV[ results of the fidelity calculation and the EE as 
well as the Schmidt gap of the system are presented. A 
discussion is provided in the last section. 


II. HAMILTONIAN 


The Hamiltonian of a ID spin-1 XXZ chain with alter¬ 
nating single-site anisotropy is given by 

N 

H = ^ J(5f + {-lyDiSn^ 

( 1 ) 

where S°‘{a = x,y,z) are spin-1 operators on the i-ih. 
site and N is the length of the spin chain. The peri¬ 
odic boundary condition is considered, and it is denoted 
that -I- 1 = 1. The parameter J denotes the antiferro¬ 
magnetic coupling, and J=1 is considered in the paper. 
The parameters A and D are the anisotropic spin-spin 
interaction and single-site anisotropy, respectively. 


III. MEASUREMENTS AND ALGORITHM 

A concept from quantum information theory, i.e., the 
ground-state fidelity, can be applied to capture the oc¬ 
currence of the QPTs. Taking a general Hamiltonian 
H[D) = Hq + DHj as an example, where Hq is the main 
part, Hj is the driving part and the quantum parame¬ 
ter D denotes its strength. If p{D) represents a density 
matrix of the system, the ground-state fidelity between 
Pq{D) and Pq[D -\- SD) can be defined as 

F{D, SD) = Tr[^ pI/^{D)po{D + SD)pI/^{D)], (2) 


where SD is a small deviation. For a pure state po = 
IV^o)(V^o|, Eq. ([|) can be rewritten as 

FiD,SD) = \{MD)\MD + SD))l (3) 

which represents the overlap of the wavefunctions at 
two adjacent quantum parameter points, and F{D, SD) 
reaches its maximum value Fmax = 1 at SD = 0. Ex¬ 
panding \ij^o{D -I- SD)) to the first order, we have 

\MD + SD)) = \md))+5dY, ( 4 ) 

y%Eo{D)-E„{D) 


where Ffno = {'ijjniD)\Hj\'ijjo{D)), and the eigenstates 
\ilJn{D)) satisfy H{D)\'tjjn{D)) = £'„|V'„(D)). Therefore, 
the fidelity susceptibility removes the artificial variable 
SD [231 US] and can be calculated by 


Xf{D) 


-2\nF{D,SD) 

{Wy ’ 


and then it yields 

[Eo{D)-E^{D)y 


(5) 

( 6 ) 


The divergence of FS can directly locate the criti¬ 
cal points. The related phase transition is exactly 
convincing [25| |26l [34] . 

The EE can be chosen as a measurement of the bipar¬ 
tite entanglement, which is defined as follows. Assuming 
\g.s.) is the ground state of the target Hamiltonian, which 
can be divided into two subsystems A and B. One con¬ 
venient choice of subsystem A is composed from the first 
site to the Lth site and the subsystem B is the rest of 
the system. The reduced density matrix of part L can be 
obtained by taking the partial trace over system N — L, 
which is given by 


PA=TrB{\g.s.){g.s.\). (7) 


Then, the bipartite EE measures the entanglement be¬ 
tween parts A and B as 


Sl = -Tr{pA log2 Pa)- (8) 

In addition, as a local order parameter, the Schmidt gap 
can also be used to describe the QPTs[35|. It is defined 
as 


G — gi — g 2 , (9) 

where gi and 52 are the first and the second largest eigen¬ 
values of the reduced density matrix pA [Eq. 0 ] , respec¬ 
tively. 

Thanks to the density-matrix renormalization- 
group(DMRG) [36113^ method, the ground state of the 
ID system can be calculated with very high accuracy. 
We implement GPU speeding up Matlab code for the 
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FIG. 1. (Color online) The fidelity susceptibility per site is 
plotted as a function of the single-site interaction for different 
system sizes with A = 1. Inset: jN for various sizes N. 

The line is a linear fit for guiding eyes. 


finite-size DMRG with double precision. The maximum 
eigenstates kept is m = 200 during the procedure of 
basis truncation, and the truncation error is smaller 
than 10“® for system sizes up to iV = 100. With such 
high performance calculation, we can precisely analyze 
the QPTs in terms of both the EE and the FS. 


IV. NUMERICAL RESULTS 

As a check, we plot the FS per site as a function of 
the single-site interaction for different system sizes with 
A = 1. As shown in the Fig. [T] peaks in the ground- 
state FS are observed, which signal precursors of the 
phase transition from Haldane phase to periodic Neel 
spin solid at Dc = 3.30, and this was confirmed by pre¬ 
vious results [33]. 

In Fig. [2];a), we plot the FS per site of XXZ chain as 
a function of single-site anisotropy D for different lattice 
sizes V=40, 60, 80. Two peaks of the ground-state FS 
are observed, indicating that the dramatic changes in the 
ground-state structure with the increase of D take place 
twice. Initially, i.e., D = 0, the system is in Neel state 
I titi ■ • • ti) when A = 2. The system becomes periodic 
Neel state jt04,0t04,---0t0 4,) when D is large 
enough. In order to characterize the intermediate phase, 
we calculate the string order parameter (SOP), whose 
definition is given by|2] 



D 

FIG. 2. (Color online) (a) The fidelity susceptibility per 
site is plotted as a function of the single-site anisotropy D for 
different system sizes with A = 2. Inset: x ™°'^for various 
sizes N. The lines are fitting lines, (b) The second derivative 
of the ground-state energy density is plotted as a function 
of the single-site interaction for different system sizes with 
A = 2. 
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FIG. 3. (Color online)The string order parameter Ox is 
plotted as a function of the single-site anisotropy for different 
system sizes with A = 2. 


Ox = — lim [5'f exp(*7r S'f)S'“]. (10) 

i<l<j 

The SOP characterizes the topological order in the Hal¬ 
dane phase[39l|40]. The results are shown in Fig[^ where 
Ox is nonzero in the intermediate phase, implying it is the 
Haldane phase. The finite-size scaling of the FS shows 
a power-law divergence at criticality, and this suggests 


that both QPTs belong to a second-order transition[25]. 
In Fig. [2](b), we also study the second derivative of the 
ground-state energy density, which exhibits divergence at 
criticalities. This confirms that the transitions should be 
of second order. 

In order to locate the QPTs precisely, the EE between 
two half subsystems is plotted as a function of D for 
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FIG. 4. (Color online)(a)The block entanglement between 
two half subsystems labeled by E is plotted as a function of D 
for different system sizes with A = 2.(b),(c) Finite-size scaling 
of Dc of the EE. The lines are the fit lines. 


different system sizes in Fig. |^a). It is observed that 
starting from small D the entanglement firstly grows as 
the D increases, meaning the quantum fluctuation en¬ 
hance the quantum nature in the classical-like system. 
After reaching a local maximum, the variation of the EE 
becomes relatively small. When the D increases further, 
the EE reaches another maximum and then decreases 
rapidly. Eurthermore, the EE is insensitive to the system 
size, as a consequence of boundary law in ID gapped 
Hamiltonian, while the maximum values of the EE grow 
with increasing Af, which is a logarithmically divergent 
correction to the boundary law[lTl 32] ■ The location of 
left peak in FigQa) moves to higher D up to a particular 
value as the system size increases, while the location of 
right peak moves to lower up to a particular value as 
the system size grows. We fit the locations of maximums 
by the formula 

D,{N) ^ D, + aN-\ ( 11 ) 

where a, b are size-independent constants and N is the 
system size. We obtain that Dd = 3.20, Oi = 14.44, 
hi = 1.438 and Dc 2 = 4.760, 02 = —14.45, 62 = 1.506; 
see Fig. ib) and (c). 

We also investigate the Schmidt gap of the reduced 
density matrix after cutting a A^-site chain into two half 
subsystems {L = N/2). The Schmidt gap labeled by G 
is plotted as a function of D for different system sizes in 
Fig. [5| It is seen that the Schmidt gap is large when 
the system is in the Neel phase. With increasing D, 
the Schmidt gap closes very rapidly when the system is 
in the Haldane phase, because S=1 Haldane phase is a 
topological phase protected by specific global symmetry 
and is characterized by a double degeneracy of the en¬ 
tanglement spectrum |43j. When D increases further, the 



FIG. 5. (Golor online)The Schmidt gap labeled by G is 
plotted as a function of D for different system sizes with A = 2. 
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FIG. 6. (Golor online)Phase diagram of a spin-1 Heisenberg 
chain with alternating anisotropy D and A. Here 0 and 
4 stand for the single-site states |S'f) with Sf = 1, 0 and -1 
respectively. t4t4t4 (t 0 4 0 t 0 4 ) denotes the Neel (periodic 
Neel) phase. The dot line D = 2A. 


Schmidt gap opens very rapidly again when system is in 
the periodic Neel phase. 

In Fig. we portray the D-X phase diagram of 
the Hamiltonian ([^, which is detected via the EE, the 
Schmidt gap and the FS. When D = 0, the quantum 
transition from Haldane spin liquid to Neel spin solid 
is continuous and belongs to a second-order QPT at 
A = 1.18[3l [ 22 ] . The critical points between the peri¬ 
odic Neel phase and the Haldane phase decrease with A 
increases. The critical points between the Neel phase and 
the Haldane phase increase as A increases. Two critical 
lines will merge each other when A increases, where the 
Haldane phase disappears and the Neel—periodic Neel 
transition becomes first order at D = 2X. 


V. DISCUSSION 

In the present paper, we have investigated the QPTs 
in the ID spin-1 XXZ chains with alternating single-site 
anisotropy by analysing the bipartite entanglement, the 
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Schmidt gap and the FS by using the DMRG technique. 
Their relation with QPTs is under discussion. It is impor¬ 
tant to note that the quantum phase transitions from the 
Neel ordering to Haldane spin liquid to periodic Neel spin 
solid can be well characterized by the FS. The finite-size 
scaling demonstrates that FS should diverge in the ther¬ 
modynamic limit at the pseudo-critical point, and the 
locations of extreme points approach quantum critical 
point accordingly. It also shows a power-law divergence 
at criticality, which indicates the QPT is of second order, 
and the critical exponent can be obtained. We compare 
the FS with the second derivative of the ground-state en¬ 
ergy, and hnd both of them exhibit similar peaks. The 


critical point can also be successfully detected by the EE 
and the Schmidt gap. To sum up, the quantum informa¬ 
tion observables are effective tools for detecting diverse 
QPTs in spin-1 models. 


ACKNOWLEDGMENTS 

This work is supported by the National Natural Sci¬ 
ence Foundation of China (NSFC) under Grants No. 
11104021, No. 11347008 and No. 11474211, as well as 
the Natural Science Foundation of Jiangsu Province of 
China under Grant No. BK20141190. 


[1] S. Sachdev, Quantum Phase Transitions (Cambridge 
University Press, Cambridge, England, 1999). 

[2] M. den Nijs and K. Rommelse, Phys. Rev. B 40, 4709 
(1989). 

[3] W. Chen, K. Hida, and B. C. Sanctuary, Phys. Rev. B 
67, 104401 (2003). 

[4] C. Degli Esposti Boschi, E. Ercolessi, F. Ortolani, and 
M. Roncaglia, Eur. Phys. J. B 35, 465 (2003). 

[5] J. Darriet and L. Regnault, Solid State Commun. 86, 409 
(1993). 

[6] W. J. L. Buyers, R. M. Morra, R. L. Armstrong, M. J. 
Hogan, P. Gerlach, and K. Hirakawa, Phys. Rev. Lett. 
56, 371 (1986). 

[7] M. Steiner, K. Kakurai, J. K. Kjems, D. Petitgrand, and 
R. Pynn, J. Appl. Phys. 61, 3953 (1987). 

[8] R. M. Morra,W. J. L.Buyers,R. L.Armstrong, 

andK.Hirakawa, Phys. Rev. B 38, 543 (1988). 

[9] Kiran Singh, Tathamay Basu, S. Chowki, N. Mahapotra, 
Kartik K. Iyer, P. L. Paulose, and E. V. Sampathku- 
maran, Phys. Rev. B 88, 094438 (2013). 

[10] A. Zheludev, J. M. Tranquada, T. Vogt, and D. J. But- 
trey, Phys. Rev. B 54, 7210 (1996). 

[11] S. Raymond, T. Yokoo, A. Zheludev, S. E. Nagler, A. 
Wildes, and J. Akimitsu, Phys. Rev. Lett. 82, 2382 
(1999). 

[12] A. Zheludev, T. Masuda, 1. Tsukada, Y. Uchiyama, K. 
Uchinokura, P. Boni, and S.-H. Lee, Phys. Rev. B 62, 
8921 (2000). 

[13] A. I. Smirnov, V. N. Glazkov, T. Kashiwagi, S. Kimura, 
M. Hagiwara, K. Kindo, A. Ya. Shapiro, and L. N. Demi- 
anets, Phys. Rev. B 77, 100401 (R) (2008). 

[14] Zhe Wang, M. Schmidt, A. K. Bera, A. T. M. N. Islam, 
B. Lake, A. Loidl, and J. Deisenhofer,Phys. Rev. B 87, 
104405 (2013). 

[15] A. K. Bera, B. Lake, A. T. M. N. Islam, B. Klemke, 
E. Faulhaber, and J. M. Law, Phys. Rev. B 87, 224423 
(2013). 

[16] Jak Chakhalian, John W. Freeland, Andrew J. Millis, 
Christos Panagopoulos, and James M. Rondinelli, Rev. 
Mod. Phys. 86, 1189 (2014). 

[17] Wen-Long You, Peter Horsch, and Andrzej M. 01es,Phys. 
Rev. B 89, 104425 (2014); Wen-Long You, Guang-Hua 
Liu, Peter Horsch, and Andrzej M. Oles, Phys. Rev. B 
90, 094413 (2014). 


[18] L. Amico, R. Fazio, A, Osterloh, and V. Vedral, Rev. 
Mod. Phys. 80, 517 (2008). 

[19] Shi-Jian Gu, Int. J. Mod. Phys. B 24, 4371(2010). 

[20] Wen-Long You, Andrzej M. Oles, and Peter Horsch, 
Phys. Rev. B 86, 094412 (2012). 

[21] Guang-Hua Liu, Wei Li, Wen-Long You, Guang-Shan 
Tian, and Gang Su, Phys. Rev. B 85, 184422 (2012). 

[22] Jie Ren and Shi-Qun Zhu, Phys. Rev. A 79, 034302 
(2009). 

[23] D. F. Abasto, A. Hamma, and P. Zanardi, Phys. Rev. A 
78, 010301(R) (2008). 

[24] P. Buonsante and A. Vezzani, Phys. Rev. Lett. 98, 
110601 (2007). 

[25] W. L. You and Y. L. Dong, Phys. Rev. B 84, 174426 

( 2011 ). 

[26] W. L. You, Y. W. Li, and S. J. Gu, Phys. Rev. E 76, 
022101 (2007). 

[27] Huan-Qiang Zhou, Jian-Hui Zhao, and Bo Li, J. Phys. 
A: Math. Theor. 41, 492002 (2008). 

[28] O. Legeza, J. Solyom, L. Tincani, and R. M. Noack, Phys. 
Rev. Lett. 99, 087203 (2007). 

[29] Jie Ren and Shi-Qun Zhu, Eur. Phys. J. D 50, 103 (2008); 
Jie Ren and Shi-Qun Zhu, Int. J. Quant. Inf. 9, 531 
( 2011 ). 

]30] J. Ren, X. Xu, L. Gu, J. Li, Phys. Rev. A 86, 
064301(2012). 

[31] Y. C. Tzeng, H. H. Hung, Y. C. Chen, and M. F. Yang, 
Phys. Rev. A, 77, 062321 (2008). 

[32] Y. C. Tzeng and M. F. Yang, Phys. Rev. A 77, 012311 
(2008). 

[33] Kazuo Hida and Wei Chen, J. Phys. Soc. Jpn 74, 2090 
(2005). 

[34] M. Cozzini, R. lonicioiu, and P. Zanardi, Phys. Rev. B 
76, 104420 (2007). 

[35] G. De. Chiara, L. Lepori, M. Lewenstein, A. Sanpera, 
Phys. Rev. Lett. 109, 237208 (2012). 

[36] S. R. White, Phys. Rev. B 48, 10345 (1993). 

[37] U. Schollwock, Rev. Mod. Phys. 77, 259 (2005). 

[38] U. Schollwock, Ann. Phys. (NY) 326, 96 (2011). 

[39] Hiroshi Ueda, Hiroki Nakano, and Koichi Kusakabe, 
Phys. Rev. B 78, 224402 (2008). 

[40] Wei Li, Andreas Weichselbaum, Jan von Delft, Phys. 
Rev. B 88, 245121 (2013). 

[41] Michael M. Wolf, Phys. Rev. Lett. 96, 010404 (2006). 



6 


[42] D. Gioev and I. Klich, Phys. Rev. Lett. 96, 100503 
(2006). 

[43] F. Pollmann. A. M. Turner, E. Berg, and M. Oshikawa, 
Phys. Rev. B 81, 064439 (2010). 



